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Abstract 

The dynamical evolution of a two-site Hubbard model is derived 
in the presence of an uniform external magnetic field for a general 
initial state. The time evolution of the half-filled (two-particle) case 
has a complex behaviour. Under certain initial conditions, the average 
number of fermions in one-particle state, in the half-filled case, is not 
necessarily periodic, even though the magnetization remains periodic. 
^ ' The results obtained may be applied to study the magnetization and 

the transition dipole moment of the organic charge-transfer salts in 
the in-phase mode. 
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Dynamics of the Two-Site Hubbard Model 



Even the simplest fermionic many-body model in Condensed Mat- 
ter Physics, is too complicated to be treated exactly. This is the case 
of the well known Hubbard model Despite of its simplicity, it 
gives good qualitative descriptions of many important phenomena in 
Condensed Matter Physics. These results come from numerical simu- 
lations (of the Monte Carlo type, for example) or from a perturbative 
treatment of one of the terms in the hamiltonian |l| . Non-perturbative 
analytical treatment of the model turns out to be feasible if one fur- 
ther simplification is made: to consider a small number of space lattice 
wites. Although this simplification seems to impose serious limita- 
tions on the applicability of the model, it still can describe interesting 
phenomena occurring in real systems, e.g., the appearance of extra 
excitation lines in the infra-red spectra of organic charge-trasfer (CT) 
salts HI, ||, §]. In the literature, this system is described by two-site 
Hubbard model[||, ||]. Usually, it has been studied energy- levels || and 
the optical properties of the CT salts in the presence of the external 
time dependent electric field Q, in the half- filled case. There are also 
experimental data of the energies of the peaks in the dimer spectrum 
for different temperature. In all previous work, the system is not 
studied in the presence of an external magnetic field. 

The time dependence of the CT salts results is a consequence of an 
external electric field. However, it is missing the dynamics of the CT 
salts under general initial conditions of the vector state that describes 
the dimer. 

In this work we present the exact time evolution of the two-site 
Hubbard model, with arbitrary band filling and coupling constant, in 
the presence of an uniform magnetic field. From an arbitrary many- 
body initial state we obtain the exact density operator of the system. 
The one-body properties of this identical particles system are obtained 
by calculating the one-particle reduced density matrix. These results 
give us a dynamical picture of the behaviour of the one-electron states 
of the system, and allow us to obtain the electron population in each 
one-particle state, the magnetization per site and the transition of 
the electric dipole moment of CT salts in the in-phase mode, in the 
presence of an external constant magnetic field. 

The well-known Hubbard model |l[ has the second-quantized hamil- 
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tonian 



ij,<r i,a i 

where a.j )CT is the destruction operator of one electron with spin a 
on Sit© % 5 t<i j are the so-called hopping integrals, U is the effective 
intra-atomic Coulomb interaction, Xb = — 3M f — , where 5 is the Lande 
factor, fis is the Bohr magneton, is the external magnetic field 
chosen in the z direction, and nj j(T = aJ CT aj )(J . The creation (aj CT ) and 
destruction (a., jCr /) operators satisfy the anti-commutation relations: 

{ a l,<x> a .?v} = s ij s ™' > 

(2) 

{a iiff ,a ii(T ,} = {af^at^,} = 0. 

We consider here the case where the hamiltonian ([!]) takes into 
account only nearest-neighbour hopping. We also have tn = £22 = Eq 
and t\2 = t2i = T. Under these conditions, a two-site Hubbard hamil- 
tonian ([l]) describes, for example, the CT salts in the in-phase mode 
[||, but in absence of vibronic couping with the internal modes of 
the monomers. For the CT salt, £"0 denotes the enrgy of the radical 
electron molecular orbital, T gives the hopping integrals of electrons 
between the molecules and U is the effective Coulomb interaction be- 
tween the two elctrons on the same molecule ||]. 

In order to study the most general state of the system for the two 
sites case, we choose a basis of eigenstates of the number operator 
n^tj, whose components can be written as 

|n 1 n 2 ;m 1 ,m 2 ) = (a^)** (a^r (a 24 ) mi HtH 1°) > ( 3 ) 

where nj, rrtj = 0,1 and i = 1,2. |0) represents the vacuum state of 
the model. 

The action of H over each of the basis states can be obtained by 
straightforward calculation, using the anticommutation relations eq. 

(D- 

The most general initial condition for the system described by a 
state vector is 
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1 

1,712, mi, m2 

(ti)\n 1 n 2 ;mim 2 ) , (4) 

ni ,712, mi ,m2=0 

where \n\n2\vn\m2) are given by eq. (|3j) and the values of the constants 
/ni,7i 2 ,mi,m 2 (0) are determined by the initial conditions. Due to the 
fact that the eigenstates of the number operator, eq. (|3|), form a 
complete basis, at any time we have, 

1 

1,712,5711,7712 

(t)\nin 2 ; mim 2 ) ■ (5) 

ri\ ,ri2,mi ,m2=0 

Inserting (^(t)) in the Schrodinger equation gives a system of cou- 
pled first order differential equations for the coefficients fni,n2,m\,rm(t) 
that can be solved analytically. Actually, before solving the cou- 
pled equations for the coefficients, it is worth to notice that, since 
[H, N CT ] = 0, where N CT = J2i n o-,i, this system of coupled equations 
brakes up in smaller set of systems. The dynamics of states with dif- 
ferent spin components are decoupled. We end up with five systems 
of first order differential equations. In appendix A, we give the ex- 
plicity time evolution of each eigenstate of the total number operator 
(N = J2i o- n <T,i)- Not all vectors belonging to basis ([|) are eigenstates 
of the hamiltonian (jl]), but since N and N CT are constants of motion, 
the eigenstates of H are labeled by their eigenvalues of N and N CT . 

With the time-dependence of the most general state of the sys- 
tem in hands, we can write the time-dependent density matrix of this 
many-body system, 

P (t) = m))(m\- (6) 

If we are interested in the effective dynamics of the one-particle sub- 
system, the one-body density matrix ||, ||| must be determined. It is 
defined in terms of the full density matrix as 

Ai,a;jy(t) =Tr{ a ] <T , ai>a p(t)}. (7) 

For the sake of simplicity, instead of working with four indices in 
the definition of A(t) we redefine them as: 

a U = &l> a 2,| = a 2 , a^f = a 3 , a 2 ,t = a 4 , (8) 
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and in an analogous way the creation operators. 

Using the new definition for the indices, and, the fact that p(t) is 
a pure density matrix, the trace in eq.(0) reduces to 

Aitf) = (*(t)|atai|*(t)> i,j = !,-■■, 4- (9) 

It is clear from the definition ([?]) that the diagonal elements of A(t) 
give the average population of the one-particle state i for the identical 
particle system described by ^(i)). From these elements, other im- 
portant quantities may be determined, such as magnetization per siteQ 
and electric dipole moment^. It is interesting to notice that, although 
the whole system is in a pure state, the effective one-body subsystem, 
depending on the initial values of /ni,n2,mi,m 2 (0)> * s m a statistical 
mixture of states. This may be easily verified by direct calculation of 
A 2 (t) and subsequent comparison with A(t) itself. In general, we get 
A 2 (i) 7^ A(t), in a clear indication of the statiscal mixture character 
of the one-particle subsystem state. 

Although we have an analitical expression for the most general 
state of the system \^f(t)), it is very cumbersome and not worth show- 
ing here (look in Appendix A for the time dependence of the coeffi- 
cients fmimwxnzit) ec l- ©)■ Instead, we choose a set of numerical 
values for the constants of the model and some initial states ^(O)), 
and present the results for some one-particle quantities. 

The zero- and four-particle sectors of the Fock space are very un- 
interesting. They have only one vector each, which evolves in time 
only by a phase. If the initial state |^(0)) lies entirely in one of these 
sectors, all the one-particle observables are constant in time. 

The one- and three-particle sectors have both four vectors, and are 
divided into two sets of two coupled vectors each. If |^(0)) lies within 
one of these sectors, the one-particle observables are time-dependent 

1 The z component of the magnatization of site i is|| 

m- = -s/iB(n u -n u ), 

where g is the Lande factor and \xb is the Bohr magneton. 

2 The electron electric dipole moment of the dimerj^]: 



where = X^=t,i n '-^- 



P = y(ni - n 2 ), 



in general, being periodic functions of time with the frequency 

Since the dynamics does not couple different Fock sub-spaces, the 
diagonal elements of the A matrix in these two Fock-subspaces are: 



A-11 = | /(1,0;0,0;0) | 2 cos 2 (— )+ | /(0,0;1,0;0) | 2 sin 2 (— ), 



(12) 



Ag = | /(1,0; 1,1; 0) | 2 + | /(l, 1; 1, 0; 0) | 2 + 

Tt Tt 
+ | /(l, 1; 0, 1; 0) | 2 cos 2 (— )+ | /(0, 1; 1, 1; 0) | 2 sin 2 (— ) 



(13) 



and 



Ml = | /(0,1;0,0;0) | 2 cos 2 (— )+ | /(0,0;0,1;0) | 2 sin 2 (— ), 



(14) 



Ag = |/(0,1;1,1;0) | 2 + |/(1,1;0,1;0) | 2 + 



Tt Tt 
+ | /(l, 1; 1, 0; 0) | 2 cos 2 (— )+ | /(l, 0; 1, 1; 0) | 2 sin 2 ( T ), 

(15) 

where A^' are diagonal elements of the one-particle reduced density 

(3) 

matrix in the Fock sub-space N = 1 and A^ the diagonal elements 
for N = 3. In the solid-state language, these two sectors correspond to 
the quarter and three-quarter band filling of the model respectively. 
Due to the particle-hole symmetry of the model, these two sectors 
have very similar characteristics. Note that when the state has three 
fermionic particles, there is interaction among them. However, we see 
from eqs. (14) and (|l5|) that the dynamics of theses states is only 



dictated by the the hopping integral T, as is compatible with particle- 
hole symmetry. 

The operator N ff is a constant of motion. To describe the dynamics 
of the one-particle sub-system it is enough to discuss, for fixed a, the 
population in one of the space sites. We have chosen to fix in eqs. 
(|l|)-(0) and eqs.(|l§-(|l5|) on site % = 1. 
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The two-particle sector has the richest structure, for it has six 
state vectors, where four of them are coupled. If we choose |^(0)) 
to lie in this sector, the one-particle observables are in general quasi- 
periodic functions of time. The time dependence of the population 
of each one-particle observable has two frequencies: W\ = |^ and 

W2 = ^ U ' 2 2h 6T ~ ■ ^ n § enera b where p and q are integers, which 

implies that the dynamics of the population of the one-particle sub- 
system is not periodic in general. For some special initial states, only 
a few of the one-particle observables turn out to be periodic. This 
happens because of subtle cancellations occurring in the calculations 
of those observables, as will be shown later on. This sector correspond 
to the half-filled Hubbard model, widely studied in solid state physics. 
In particular, this is the band-filling used to describe the charge-trasfer 
organic salts j|, ||] . 

In order to ilustrate the qualitative discussion given above, we 
present some graphs for one-particle observables obtained from the 
analytical expressions. We need to point out that the external con- 
stant magnetic field only appears in the phases of the time evolution 
of states @ of the Fock sub-spaces N = 1, N = 3 and for states of 
N = 2 when both fermions have the same spin component. But, in 
the one-particle reduced density matrix, it gives no contribution. Let 
us consider the graphs when the model constants are chosen as: 



T = 1 
U = 0.5, 

X B = 0. (16) 

We have taken h = 1. In doing this, T/Ti turns out to be an energy 
scale for the system. We first chose U comparable to the hopping in- 
tegral T to see how the two terms in the Hamiltonian combine when 
none of them is dominating. This competition reveals itself in the 
time-dependence of the one-particle observables by the absence of a 
characteristic frequency in their oscillations. In figure [l] we show the 
average spin up occupation = A-i^ii = A^t), and the average 
electric dipole moment of the dimmer d(t) = 2i(t) + A2 i 21 (^) — 
(A ltjlT (i) + A u , u (t)) = A 22 (t) + A 44 (i) - (A 33 (i) + A n (t)) for the 
initial state |^i(0)) = |11)|00). It is easy to identify both functions 
as non-periodic by simply looking at the figures. The average magne- 
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tization per site of this state is zero, since its initial value is zero and 
it is a constant of motion. 

If we now choose a new initial state |* 2 (0)> = ^(|10)|01) - 
|01)|10) + [11) |00) ) , we obtain two non-periodic observables, namely 
n\i(t) and di(t), but the average magnetization of site 1, m\{t) = 
A-lt,lt(*) ~~ ^-11,11 (*) j turns out to be periodic ( see Figure §). This oc- 
curs because of a cancellation of the contributions from 1 10) 1 01) and 
|01) 1 10) to the average magnetization per site, due to the fact that 
both states enter the combination with coefficients of same modulus 
but opposite sign. 

In order to ilustrate how the relation between the hopping integral 
T and the Coulomb interaction U may influence the observables of 
the system, we choose a new set of model constants, and repeat the 
calculations already performed. In the literature ||] it is stated that for 
the CT salts U = 5.45T, therefore, we choose the new set of constants 
as: 

T = 1 
U = 5.45, 

Xb = 0. (17) 

Now the Coulomb interaction U is four times greater than the hop- 
ping integral T. In figure ^ we observe that this noticeable difference 
between them causes the appearance of a modulation of frequencies. 
The higher frequency oscilations due to the Coulomb integral are su- 
perimposed to the hopping integral oscilations. 

In summary, we studied the exact dynamics of two-site Hubbard 
model when the identical particles system is described by a vector 
state. We considered the case where t%2 = £21 = T, that has been 
used as a model to explain the CT salts in the in-phase mode Q, 
without the interaction of the radical electrons with the internal vi- 
bronic modes of the molecule. 

The particle-hole simmetry of the Hubbard model is well realized 

(3) 

by the time evolution of A- , i = 1,2. Even though for the three- 
particle Fock sub-space the Coulomb interaction (|l]) still works, the 
dynamics of the population of one-particle states is governed by the 
hopping integral T, as happens for the Fock sub-space N = 1. 
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In the half-filled case, depending on the initial conditions, we can 
get non-periodic behaviour for some one-particle properties, while the 
magnetization per site remains periodic. 

Finally, we have shown as the one-particle observables as magne- 
tization and site occupation oscillate in time. It is well known that, 
when these quantities are measured in real systems as magnetic tran- 
sition metals for example, the measurements furnish constant frac- 
tional values. From the theory of intinerant magnetism we have that 
the measured fractional charge and magnetization of such metals is 
a consequence of a time average performed by the measuring appa- 
ratus. We obtain fractional values by taking the time average of the 
time-dependentt observables we have calculated. This time average is 
defined as 

{w(t))= lim I [ w(t)dt, (18) 
T-*oo 1 Jq 

For the half- filled case, we get from eq.(|l8|) and what was discussed 
previously, that the average value for the magnetization per site is a 
fractional constant, but the eletric dipole transition is not, since it 
is not a periodic function in time. We get a time dependent electric 
polarizability in the absence of an external electric field. It is an open 
question if this effect can actually be measured. 
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Appendix A: Time Evolution of the Eigenstates 
of nj )Cr Operator 

In this apendix we present the explicit time dependence of the 
most general state of the two-site Hubbard model. It is obtained by 
inserting 

1 

im 2 ni?i2 (i)|mim 2 )|nin 2 ) (19) 

mi,m2,ni,ri2=0 

in the Schrodinger equation, and then projecting out this equation over 
each one of the basis states This will, in principle, give a set of 
sixteen coupled first order differential equations for the f mi m2n 1 n 2 (t) ■ 
Due to some simetries of the Hamiltonian, these sixteen equations de- 
couple, forming a set of five independent systems with variable number 
of equations. From the hamiltonian ([!]), we have the foiling constants 
of motion: 



[H,N CT ] = 0, 
[H,N] = 0. (20) 



Eq.(20) leads to total charge and spin conservation of the hamiltonian, 
which defines five independent sectors of the Hilbert space, labeled by 
the total charge of the states. These sectors may be further divided 
in subsectors of states with the same total spin (or magnetization). 
Therefore we may determine the dynamics of each subsector indepen- 
dently. The explicit time dependence of each coefficient f mim2 n 1 n 2 (t) 
is given below. 



Fock sub-space N=0: 

/oooo(i) = /oooo(0) 

Fock sub-space N=l: 

Subsector a =[: 

£ ri /an /*12*21,n .\t\ , , . . ,£l2*21 ,ni _i (e ~ A -B ) 

/loooW = L/iooo(0)cos(— ^— t) - i— /ooio(0)szn(— — *)Je fi 

£ u\ rj m\ /*12*21,s . \t\ . , . . ^12*21 _i (e ~ A -g) 

/ooio(i) = L/ooio(0)cos(— — t) - i — /iooo(0>m( t)\e * 

a t a 
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Subsector a =]: 

/oioo(t) = [/oioo(0)c OS (^t) - iM/oooi (OjamC^OJe-*^* (24) 

ft t a 

/oooi(t) = [/oooi(0) COS (^t) - ;M /oioo (o) s ;n(^t)]e-*T' (25) 

Fock sub-space N=2: 

Subsector cr = 0: 

/noo(t) = *^gS*n(^)(/oiio(0) - /iooi(0)) + ^[/Wf <) " ^sm(^)] x 

x(r/oon(0) + r*/ 11 oo(0))e- t ^ t - ^e' 2 ^' (rf O on(0) - r*/noo(0)) 

(26) 

/oon(t) = »^«n(^)(/oiio(0) - /iooi(0)) + -j-[/?cos A) - if/sm A)] x 

■ (2e+U) . I r,-(£+U). 

x(r/oon(0) + r*/iioo(0))e- t ^ t + -e" 2 *— *( r / oii(0) - r*/noo(0)) 

zr 

(27) 

/iooi(t) = ^(/ono(0) + /iooi(0))e- 2 ^* - {^[r/oon(0) + r*/iioo(0) + 
+|(/oiio(0)-/iooi(0)]sin(^)- 

-^(/ono(0) - /iooi(0))co S (^)}e- i T* (28) 

/ono(t) = ^(/ono(0) + /iooi(0))e- 2i t* + {^[r/oon(0) + r*/noo(0) + 
+|(/ono(0)-/iooi(0)]sm(^)- 

-^(/ono(0) - W0))cos(^)}e-*T* (29) 
Subsector cr = — 1: 

/ioio(i) = e- 2l ^Vioio(0) (30) 
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Subsector a = 1: 

/oioi(t) = e- 2i ^^*/oioi(0) 



Fock sub-space N=3: 

Subsector a =[: 

Irl Irl Irl .(3e+2t/-A s ) 

/ino(i) = [i^r/ion(0)sm(y*) + /mo(0)c s( V^)] e "* R 

r* ft ft 

Irl Irl Irl .(-m+2u-\ b ), 

/ioii(t) = * U /nio(0)am(Vt) + /ion(0)c OS (V i )] e " 1 1 

r ft ft 



Subsector cr =|: 

Irl Irl Irl .(3n+2u+\ B ) 
/noi(t) = [^/oin(0)sm(^t) + /noi(0)co S (^t)]e- t R 

t| „ . Jrl N „ /r , ,W\ , n , (3e+2t/+A B ) 



/oin(i) = [* — /noi(0)sm(— t) + fom{0)cos(—t)]e 

T ft ft 

Fock sub-space N=4: 

/mi(t) = e- 4 ^ t /iiii(0) 

Here, fm 1 m 2 n 1 n 2 (Q) is the coefficient of the |mra2)|nin2) state at i = 0, 
r = t 12 = t* 21 and /? = ^[/ 2 + 4|r| 2 . 
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Figure 1: Average occupation of site 1 (a) and electric dipole moment (b) 
for the first set of chosen constants flip] ) and initial state |\l/i(0)) = |11)|00). 



Figure 2: Average occupation (a) and magnetization (b) of site 1 and elec- 
tric dipole moment (c) for the first set of chosen constants and initial state 
|^2(0)). Observe the non periodicity of (a) and (c), in contrast to (b), which 
is clearly periodic. 



Figure 3: Average occupation of site 1(a) and electric dipole moment (b) for 
the second set of chosen constants ([T7|) and initial state |^ 2 (0)). Contrary to 
what was seen in figs. [3] and ||, these observables present two clearly distinct 
characteristic frequencies superimposed, due to the discrepancy between the 
hopping integral and the Coulomb interaction. 
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